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abstract 

We develop perturbative QCD factorization theorems for the exclusive 
semileptonic heavy baryon decay A& — > plu, in which the relevant hadronic 
form factor is expressed as the convolution of a hard subamplitude with 
the Ab baryon wave function and the proton wave function. The specific 
evolution scale for the proton wave function, determined from the best fit to 
the data of the proton form factor, is adopted. The Sudakov resummation 
for a heavy-light system (the A& baryon) and the quark-level decay diagrams 
with at least one hard gluon attaching the b quark, which were ignored in the 
literature, are incorporated. It turns out that these additional ingredients 
are important: the neglect of the former and the latter reduces the results of 
the form factor by factors 1.5 and 3, respectively. We present the predictions 
of the form factor and of the proton energy spectrum for the various choices 
of the parameter involved in the A& baryon wave function. 
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I. INTRODUCTION 



Recently, we have applied perturbative QCD (PQCD) factorization the- 
orems to various B meson decays, and progresses have been made [1-8]. For 
exclusive processes |], ||, the form factors involved in decay spectra and de- 
cay rates are factorized into the convolution of a hard subamplitude with 
meson wave functions. The former is process-dependent and calculable in 
perturbation theory at the parton level. The latter, absorbing long-distance 
dynamics of B meson decays, are universal (process-independent), and must 
be extracted from experimental data or derived by using nonperturbative 
methods. For example, we have determined the B meson wave function from 
the best fit to the experimental data of the branching ratio B(B —>■ K*j) [§]. 
The resummation ||, [T(J of double logarithms from the overlap of collinear 
and soft enhancements in radiative corrections to meson wave functions plays 
an important role in the analyses. The resultant Sudakov factor suppresses 
the long-distance contributions, and improves the applicability of PQCD 



around the few GeV scale Wl 



The Wilson evolution in effective field theory has been incorporated into 
the above PQCD formalism, leading to the three-scale factorization theorem 
for nonleptonic decays 0. As the evolution scale runs to below the b quark 
mass Mb and the c quark mass M c , the constructive and destructive interfer- 
ences between the external and internal W-emission contributions involved 
in bottom and charm decays, respectively, appear. Nonfactorizable and non- 
spectator contributions can be evaluated systemmatically, which make possi- 
ble the simultaneous explanation of the ratios R = B(B —>■ J/ipK*)/B{B — ► 
J/ipK) and R L = B(B -> J/^K* L )/B(B -> J/ipK*) in the B -> J/ipK^ de- 
cays fl2"f . With the further inclusion of the all-order soft gluon exchanges, the 
mechanism of the opposite signs of nonfactorizable contributions in bottom 
and charm decays has been understood [|J. 

Similarly, inclusive decay spectra and decay rates of the B meson are 
factorized into the convolution of a hard b quark decay subamplitude with 
a B meson distribution function and several jet functions for light energetic 
final states || [J. The distribution function, being the outcome of the re- 
summation of nonperturbative power corrections to the processes, has been 
extracted from the photon energy spectrum of the decay B — > J7| . The 
jet functions collect the double logarithms appearing at the end points of 
decay spectra, whose resummation gives the Sudakov factor. The three-scale 
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factorization theorem also applies to inclusive nonleptonic decays. Using the 
above formalism, we have observed that the single-charm mode b — > cud is 
enhanced more than the double-charm mode b — > ccs is. The semileptonic 
branching ratio -Bsl = BiB — > Xlv) is then reduced without increasing the 
charm yield n c per B decay, and the large -Bsl controversy in the conven- 
tional approach based on heavy quark effective theory is resolved ||. Our 
predictions for the absolute lifetimes of the B meson and of the A& baryon 
are also consistent with the data, such that the lifetime ratio r{Kb)/r{Bd) 
can be explained ||. 

With the above successes, it is natural to extend the PQCD formalism to 
more complicated heavy baryon decays. In this paper we shall start with the 
simplest exclusive mode A& — > plv, taking into account the Sudakov resum- 
mation for a heavy-light system (the A& baryon in this case), investigating 
the applicability of PQCD to heavy baryon decays, and understanding the 
sensitivity of predictions to the variation of the A& baryon wave function, 
which remains unknown. It is expected that a PQCD analysis of heavy 
baryon decays is not as reliable as that of heavy meson decays, since partons 
in the former case are softer and Sudakov suppression is weaker. We shall 
show that PQCD for the decay A& — > plv is reliable only at the high end of 
the proton energy. Therefore, our approach may be appropriate for exclusive 
nonleptonic heavy baryon decays, which will be studied in the future. 

The first attempt to apply the PQCD factorization theorem to the A& — > 
plv decay was made by Loinaz and Akhoury ||13|| . However, they did not 
consider the Sudakov resummation for the A& baryon, and their evaluation 
of the hard subamplitude was not complete: the quark-level decay diagrams 
with at least one exchanged gluon attaching the b quark were ignored. Our 
analysis shows that these two ingredients are in fact important. The neglect 
of the former and the latter reduces the results of the relevant form factor 
by factors 1.5 and 3, respectively. We shall adopt the King-Sachrajda (KS) 
model [TJj for the proton wave function, along with the specific choice of its 
evolution scale, which were determined from the best fit to the experimental 
data of the proton form factor \TE, IE]. Furthermore, the full expression of 
the Sudakov factor with the accuracy up to next-to-leading logarithms [0 
will be inserted. 

In Sect. II we briefly explain how to factorize the perturbative and non- 
perturbative contributions to the decay A& — > plv into the hard subamplitude 
and the baryon wave functions, respectively. The large logarithms contained 
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in these convolution factors are organized in Sect. III. The factorization for- 
mula for the relevant form factor is presented in Sect. IV. Numerical analyses 
are performed in Sect. V, where the behaviors of the form factor and of the 
proton energy spectrum are exhibited. Section VI is the conclusion. 



II. FACTORIZATION THEOREMS 

The amplitude for the semileptonic decay A& — > plv is written as 

M = ^K 6 AUW-75h], (1) 

where Gp is the Fermi coupling constant, and V ub the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element. All QCD dynamics is contained in the 
hadronic matrix element, 

M, = (P(p')\u^(l- l5 )b\A b (p)) , (2) 

where p (p') is the momentum of the A& baryon (proton). The general struc- 
ture of a matrix element of charged weak currents between baryon states is 
expressed as 

(B',p',s'\f\B,p,s) = B'(p', S ')[f 1 (q 2 )r-if2(q 2 )cT^qv 

+f 3 (q 2 )q"]B(p,s), 
(B',p', S '\j»\B,p, S ) = B'{p\s')[ gi {q 2 )^ lb -ig 2 {q 2 )a^q v 

+93(q 2 )l^}B(p,s) , (3) 

where is the vector current, jg the axial vector current, B(p, s) the spinor 
of the initial-state baryon \B,p, s) with p and s the momentum and spin, 
respectively. The spinor B'(p', s') is associated with the final-state baryon 
{B',p', s'\. The form factors fi(q 2 ) and gi(q 2 ) depend only on the momentum 
transfer q 2 = (p — p') 2 . 

In the case of massless leptons with 

9^(1-75)^ = 0, (4) 

it is easy to observe that only the four form factors f\, f 2 , <?i, and g 2 con- 
tribute. We rewrite the hadronic matrix element as 

M, = P(p')[L 1 (q 2 h,(l- l5 )+R 1 (q 2 ) 1 ,(l+ l5 ) 

+L 2 (q 2 )ia^(l - lb )q» + R 2 (q 2 )ia^(l + ^)q u ]A b (p) , (5) 
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with the definitions 

Li = ^(h + gi), Rx = \{h-9i), 

L2 = -\(f2+92), R2 = -\(f2-92). (6) 

Since the form factor L\ dominates for unpolarized baryon decays, we 
shall evaluate only L\ in the present work. 

The hadronic matrix element involves both nonperturbative and pertur- 
bative dynamics. To evaluate it, we employ factorization theorems, in which 
the former is separated from the latter, and absorbed into universal baryon 
wave functions. After performing the factorization, we compute the pertur- 
bative part, i.e., the hard subamplitude, order by order reliably, and convo- 
lute it with the baryon wave functions. Certainly, this factorization picture 
should fail at some power of M\ b being the A b baryon mass, where 

nonfactorizable soft gluon exchanges do not cancel exactly. We have derived 
the leading-power factorization theorem for heavy meson decays in [|lj]. This 
formalism can be extended to heavy baryon decays straightforwardly, whose 
basic ideas will be briefly reviewed below. 

The lowest-order diagrams for the decay amplitude at the quark level 
are exhibited in Fig. 1, where the b quark is represented by a double line, 
and the symbol x denotes the electroweak vertex, from which the lepton 
pair emerges. The two gluons are hard, when the proton recoils fast, which 
give the necessary momentum transfer to make the outgoing quarks move 
collinearly and form the proton. The diagrams Figs. 1(g)- l(n) were not 
considered in |L3]|, because they were assumed to be suppressed by l/M\ b . 
We shall demonstrate that these diagrams in fact give contributions of the 
same order as those from Figs. l(a)-l(f). 

We then consider radiative corrections to Fig. 1, as shown in Fig. 2, where 
the bubble represents the two hard gluons, and explain how to absorb them 
into the hard subamplitude and the baryon wave functions. There are three 
leading momentum regions for radiative corrections, from which important 
contributions arise. These regions are collinear, when the loop momentum is 
parallel to an energetic light quark, soft, when the loop momentum is much 
smaller than M^ bl and hard, when the loop momentum is of order M^ b . We 
associate small transverse momenta with the valence quarks, which serve 
as an infrared cutoff. The collinear and soft contributions may overlap to give 
double logarithms. While the hard contribution gives only single logarithms. 
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Because of the inclusion of the transverse degrees of freedom, the fac- 
torization should be constructed in the impact parameter b space, with b 
conjugate to kx |I]> 0]- in the axial gauge n ■ A = 0, n being a gauge vector 
and A the gauge field, the two-particle reducible corrections on the proton 
side, like Figs. 2(a) and 2(b), have the double logarithms ln 2 (M^ b 2 ) in the 
fast recoil region and the single soft logarithms ln(fe/i), \x being the renormal- 
ization scale. Naturally, these corrections are absorbed into the proton wave 
function. The two-particle irreducible corrections, with the gluon attaching 
a quark in the baryon and a quark in the proton, like Figs. 2(c) and 2(d), 
contain only the single logarithms ln(6/i), whose effects are less important, 
and will be neglected in the following analysis. 

On the A;, baryon side, Fig. 2(e), giving the self-energy correction to the 
massive b quark, produces only soft divergences. If the light valence quarks 
move slowly, collinear divergences in Figs. 2(f) and 2(g) will not be pinched 
and these diagrams also give only soft divergences. However, there is 
probability, though small, of finding the light quarks in the A& baryon with 
longitudinal momenta of order M\ b . Therefore, Figs. 2(f) and 2(g) may still 
contribute collinear divergences. In conclusion, reducible corrections on the 
Ab baryon side are dominated by soft dynamics, but contain weak double log- 
arithms with the collinear ones suppressed. These corrections are absorbed 
into the A& baryon wave function. The remaining important contributions 
from Fig. 2, with the collinear and soft enhancements subtracted, are then 
dominated by short-distance dynamics characterized by the single logarithms 
ln(MA b //i), and absorbed into the hard subamplitude. 

The kinematics is as follows. The momentum of the A& baryon at rest 
is p = (p + ,p~,0) with p + = p~~ = Ma 6 /v^2. The momenta of the valence 
quarks in the A& baryon are 

fci = (p + , xtp~, k 1T ) , k 2 = (0, x 2 p~, k 2T ) , k 3 = (0, x 3 p~, k 3T ) , (7) 

where k\ is associated with the b quark, and Xi are the momentum fractions. 
The momentum of the proton, which recoils in the plus direction, is chosen 
as p' = {p' + 1 0, 0) with p' + = pp + . The parameter p, < p < 1, defined by 



is related to the energy fraction of the proton. The invariant mass of the 



P 




(8) 



lepton-neutrino pair is then given by mf 



i — 



(1 — p)M\ . The momenta of the 
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valence quarks in the proton are 

^1 = ( X lP + 5 0' ^1t) ) ^2 = ( X 2P + ? 0' ^2r) > ^3 = ( X 3P + ' 0; k 3r ) . (9) 

All the transverse momenta are assumed to be much smaller than M\ b , since 
the valence quarks are near or on the mass shell. Otherwise, the quark lines 
should be absorbed into the hard subamplitude. In this sense the transverse 
momenta play the role of a factorization scale, above which PQCD is reli- 
able, and below which QCD dynamics is parametrized into the baryon wave 
functions. For a similar reason, the fraction X\ should be close to unity, as 
reflected by the behavior of the baryon wave function in Sect. IV, such 
that k\ is roughly equal to M%. 

According to the above explanation, the factorization formula for the 
form factor L\ is expressed as 

L\ = [dx'][dx\ ^[(fbl^fpa'^yix^b^p',!!) 

x# a W^(4 x i} h, M Ab ,^ Abaf3l (x u b i7 p, //) , (10) 

with the symbols 

idx] = n d^s (1 - e x^j , [d%] = n ^ , (ii) 

where the ^-function is the consequence of the momentum conservation. The 
symbol with a prime, [dx'j, is defined similarly. Note that there are no 
the variables b\. Since we shall drop the transverse momenta carried by 
the internal quarks in the hard subamplitude, the transverse separations 
among the valence quarks of the A& baryon are equal to those of the proton. 
It has been shown that this approximation, simplifying the analysis, gives 
reasonable predictions for the proton form factor [|15| . 



The proton distribution amplitude \I/p is defined, in the transverse mo- 
mentum space, by 



1=1 

MCU^P^V^p'^) + (p"7 5 C) a/3 P 7 A(A;:,p',/i) 



8x/2A c 

- (a^CUiYl.P^nklp'^)} , (12) 
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where iV c = 3 is the color number, \P) the proton state, u and d the quark 
fields, a, b and c the color indices, and a, (5 and 7 the spinor indices. The 
second form shows the explicit Dirac matrix structure WJ\ with the normal- 
ization constant fp, the proton spinor P, the charge conjugation matrix C, 
and a^ v = [7^,7^/2. Using the permutation symmetry and the constraint 
that the total isospin of the three quarks is equal to 1/2, the three functions 
V, A, and T are related to a single function \l/ through [[HJ 

V(k[,k 2 ,k 3 ,p',fi) = - [&(k' 2 ,k[,k' 3 ,p,n) + ^(k[,k 2 ,k 3 ,p',ii)} , 
A(k[,k 2 ,k' 3 ,p',ii) = - [^(k 2 ,k[,k' 3 ,p',fi) - ^(k[,k 2 ,k' 3 ,p', fi)} , 
T(fc;,^,^,^ / i) = ^[vI/(A;;,^,^,^ / i) + vI/(^,^,A;;,^ / i)] . (13) 

The structure of the A b baryon distribution amplitude \E , a 6 is simplified 
in the heavy quark limit. Under the assumptions that the spin and orbital 
degrees of freedom of the light quark system are decoupled, and the baryon 
is in the ground state (s-wave), ^a,, is expressed as [[3~3|1 

V^afr = ^^[{i> + M Ab ) l5 C}^A ha ^{k 1 ,k 2 ,k 3 ,p,^ , (14) 

where f\, is the normalization constant, and is the Af, baryon spinor. 
There is also only a single wave function $ associated with the A;, baryon. 



III. SUDAKOV RESUMMATION 

The large logarithms, appearing in the wave functions and the hard sub- 
amplitude, should be organized by the resummation technique and the 
renormalization-group (RG) method. The RG summation of the single log- 
arithms ln(MA 6 //i) and ln(6/z) leads to the evolution from the characteristic 
scale of the hard subamplitude to the factorization scale 1/6. The resumma- 
tion of the double logarithms ln 2 (MA 6 6) exhibits Sudakov suppression in the 
large-6 (long-distance) region, such that PQCD analysis becomes relatively 
reliable at the energy scale of M\ b . The standard derivation of a Sudakov 
factor is referred to |9|, [TCf . 
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The result of the Sudakov resummation for the proton wave function \l/p 
is given by 



exp 



J2s(cw } x' l p' + ) 



1=1 

xif>(x' i , cw) 



w ft 



7g(«s(/i)) 



(15) 

-a s /n. The exponent s, 



with the quark anomalous dimension r y q (a s 
corresponding to the double-logarithm evolution, is written as JlO 

s(w,Q) 



Q dp 

w p 



In 



A(a s (p)) + B(a s (p)) 



(16) 



where the anomalous dimensions A to two loops and B to one loop are 



.4 
B 
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2 



7T 

y 



10 



27 



n/ + -/3 ln — 



3 7E 



(17) 



n/ = 4 for Ab baryon decays being the flavor number, and 'Je the Euler 
constant. The two- loop running coupling constant, 



aM = 1 _ A lnln(/i 2 /A 2 

/3 ln(/i 2 /A 2 ) /3 3 ln 2 (/x 2 /A 2 ) 



with the coefficients 



A) 



33 - 2n 



f 



0i 



153 - 19n 



/ 



(18) 



(19) 



12 ' ' 24 

and the QCD scale A = Aqcd, will be substituted into Eq. (|16D . Compared 
to |T5| , our expression of the Sudakov factor with the accuracy up to single 
logarithms is more complete, since we have used the two-loop running cou- 
pling constant in Eq. (Jl8|). The second exponent in Eq. (Jl5|), i.e., the integral 
containing 7 g , corresponds to the single-logarithm RG evolution. The func- 
tion if), obtained by factoring the p' + dependence out of coincides with 
the conventional parton model. 
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The infrared cutoff cw is the inverse of a typical transverse distance among 
the three valence quarks as stated before. The recent PQCD analysis of the 
proton form factor shows that a reasonable choice of w is [|TE| 

w = min (—, — , —\ , (20) 



bi b 2 b, 



with 6 3 = |b>i — t> 2 1 - The best fit to the experimental data of the proton 
form factor determines the parameter c = 1.14 |T(|. The introduction of 



this parameter c is natural from the viewpoint of the resummation, since 
the scale cw, with c of order unity, is as appropriate as w R ITUJ. The 



variation of c, with cw being the factorization scale, represents the different 
partitions of radiative corrections into the perturbative Sudakov factor and 
the nonperturbative wave function ip. As all 1/fej are much larger than A, 
the Sudakov form factor does not give suppression. As one of these scales 
gets close to A, the Sudakov factor tends to zero and suppresses this region. 
It is easy to observe that choosing the infrared cutoff in Eq. ( ^0|) suppresses 
all the infrared divergences. 

The Sudakov resummation for a heavy-light system, such as the B meson, 
has been developed in [|I]]. The resummation for the A& baryon wave function 
is similar, and the result is 



, , , . ^ t -no fdfi 

®{Xi,bi,p,fx) = exp 

1=2 



s(cw, xip~) - 3 / — 7,(a s (/2)) 

7 r> J CW j2 



X0(Xj, cw) . (21) 

Note the lack of the exponent s for the b quark. Since the b quark is heavy, 
radiative corrections associated with it do not produce double logarithms. 
While the coefficient of the anomalous dimension 7 g is still 3, because the 
RG evolution is determined by the ultraviolet structure of loop corrections, 
which are not affected by the heavy quark mass. The parameter c is also set 
to 1.14 for convenience. 

Working out the contraction of ^? p a ipi^H a 13 7 af3l ^A b a/3-r i n t ne momen- 
tum space, whose explicit expression from each quark-level decay diagram in 
Fig. 1 is listed in Table I, we extract the hard subamplitude H . To simplify 
the formalism, we have applied the approximation Mf, ~ M\ b , and neglected 
the transverse momentum dependence of the internal quark propagators as 
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mentioned before. The RG evolution of H in the b space is written as 



H{x' i ,x i ,b i ,M Ab ,p) = exp 



-3E 



1=1 



*j dfj, 



xH(x' i ,x i ,bi,M Ab ,t l ,t 2 ) 



(22) 



where the hard scales ti will be specified below. The two arguments t\ and t 2 
of H imply that each running coupling constant a s is evaluated at the mass 
scale of the corresponding hard gluon. 



IV. FACTORIZATION FORMULA 

Employing a series of permutations of the valence quark kinematic vari- 
ables, the summation of the expressions in Table I reduces to two terms. 
Combining Eqs. fllED , (|2l|) , and (p2|) , we obtain the factorization formula, 

2tt r 1 r 1 / A f 2n 

Li = — [dx')[dx] / b 1 db 1 b 2 db 2 / d9f P (cw)f Ab 
Z/ Jo Jo Jo 

2 

x H ji X i> X h K M^^ji^jixl, x h cw) 

xexp[-S(x' i ,x i ,cw,M Ab ,t j i)] , (23) 

where the variable 9 is the angle between hi and b 2 . The expressions of Hj 
are 



H x = a s (t n )a s (t 12 )K (J(l - x x )(l - x , l ) P M Ab b 1 )K (Jx 2 x > 2 pM Ab b 2 



H 2 = a s (t 21 )a s (t 22 )K (^x 1 x[pM Ab b 1 )K (^x 2 x 2P M Ab b 2 ) , (24) 

with Kq the modified Bessel function of order zero. The functions which 
group together the products of the initial and final wave functions, are, in 
terms of the notations, 

0123 = 0(2:i, X 2 , £3) , ^123 = ^{x[,x' 2 , x' 3 , Cw) , (25) 

given by 

^ 2(^3120132 + 2Ti 32 0l32 + ^2130123 + 2T 123 (f) 123 ) 

1 (l-x l )(l-x' 1 ) 
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24(^3120132 + 2T 132 0132 + ^2130123 + 2Ti 23 0123) 



2T: 



312^312 



+ 



-^1320312 



(26) 



+ 



2x 2 (p~ l)T 3 i20 3 12 



(1 - x 3 )(l - x' x ) (1 - x 3 )(l - x 2 ) (1 - x[)p[(l - x 2 )x[p + x 2 ] 

| [(1 - X' 2 - X 3 Xi)p - Xi]^i 32 03i2 

(1 - x' 3 )(l - xi)p{x' 2 p - 1) 
[(1 - 2xi + x 2 - x' 2 p)p + aa]T 3 i203i2 
(l-x' 3 )(l-x 1 )p(x' 2 p-l) 

[(l-2x 2 + X'x - X[p)p + X 2 ]^1320312 



+2 



(l-x 3 )(l-x 2 )p(x[p-l) 

[(1 -X 2 - X[ - (1 - X 2 )(l - X' 3 )p)p - X 2 ]T 312 (j) 312 

(l-x' 3 )(l-x 2 )p(x' lP -l) 



(27) 



The functions T are related to ip in a similar way to Eq. (|13"D. The complete 
Sudakov exponent S is written as 



S 



■) +3 / J —lq(a s (p)) 
J cw U 



1=1 



'■'ii dp 

I cw p 

"'2 (/// 
/ cu> p 



1=2 

+ J] s(cw, x' lP ' + ) + 3 T 2 (a.(/Z)) 

, i Jcto LI 



(28) 



with the hard scales 



tu = max 
t 2 \ = max 



^(1 -a;i)(l-xi)pM Ab , l/6i, 



^Jx^x[pM Ab , l/bi,cw 
tu = h2 = max ^Jx 2 x' 2 pM Ab , l/b 2 , cw 



(29) 



The first scales in the brackets are associated with the longitudinal momenta 
of the exchanged gluons and the second ones with the transverse momenta. 
All tji should be larger than the factorization scale cw by definition as shown 
above. 



For the wave function ip, we adopt the KS model |L4|], which has been 
found to give predictions for the proton form factor in better agreement with 



the data compared to the Chernyak-Zhitnitsky model JT7J. The KS model is 
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decomposed as [19|, |20| 



ip{x h w) = ipasix^^Nj 

3=0 



a s {w) 
a s (/Jo) 



(30) 



with /xq ~ 1 GeV and ipas = I2OX1X2X3. The constants Nj, aj and bj, 
and the Appel polynomials Aj are listed in Table II. The evolution of the 
normalization constant fp is given by JT7 



fp{w) = fp(fio) 



l/(6/3o) 



(31) 



with fp(fi ) = (5.2 ± 0.3) x 10~ 3 GeV 2 . 

We do not consider the evolution of the wave function <j) in the scale cw , 
because it is still not known yet. Similarly, the running of the normalization 
constant f\ b is also neglected, which is assumed to be f\ b = /p(yUo) f° r 
simplicity flI5|. It is convenient to use the new dimensionless variables, 



c 



£2 



r] = x 2 + x 3 



(32) 



X 2 + X3 

when presenting the model for the baryon wave function 0. We choose 



0(C,r / ) = iVr / 2 C(l-r ? )(l-C)exp 



2^(1 - 77) 2/3 2 77C(l-C) 



(33) 



for 12 3 and ^132, where the parameter (3 can be fixed by experimental data 
of A;, baryon decays, and m q is the mass of the light degrees of freedom. 
Accordingly, we have 



<f>(Cv) 



Nr) 2 ({l-r])(l-0 
f 2 

l Ab 



x exp 



Ml m 2 (l - r] + rj() 



2/3^(1 - C) 2/3 2 rK(l - v) 
for 03i2- The normalization constant N is determined by 



(34) 



1 rl 

d( / r)dr)<f)((,r)) 
Jo 



(35) 
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In the next section we shall investigate how the predictions for the form 
factor L\ vary with the parameter (3. 

V. NUMERICAL RESULTS 

We are now ready to compute the form factor Li in Eq. (^3p numerically, 
adopting M Ab = 5.621 GeV, A = 0.2 GeV and m q = 0.8 GeV. We consider 
two typical values of the parameter (3, 1.0 and 2.0, and the corresponding 
dependences of L\ on the cutoff of b% and 62, b\ c = &2c = b c , for the proton 
energy fraction p = 0.8, 0.9, and 1.0, are displayed in Figs. 3 and 4, respec- 
tively. It is observed that the rise of all the curves indeed saturates at about 
b c = 0.85/A, indicating that the large-6 region does not contribute because 
of Sudakov suppression. Figure 3 shows that for p = 0.8, approximately 47% 
of the contribution to L\ comes from the region with a s (l/b c )/7t < 0.5, or 
equivalently, b c < 0.61A. For p = 0.9, 50% of the contribution is accumulated 
in this region. As p = 1.0, the small-6 contribution reaches 52%. 

For the larger (3 = 2.0, which corresponds to a broader A;, baryon wave 
function 0, the percentages of perturbative contributions for p = 0.8, 0.9, and 
1.0 increase a bit up to 49%, 52%, and 55%, respectively, as shown in Fig. 4. 
Since the long-distance contribution is less enhanced by 0, the magnitudes 
of are reduced to about 63% of those for (3 = 1.0. If regarding the region 
with a s /n < 0.5 as being perturbative |IJ, the PQCD analysis of the decay 
A;, — > pW for p > 0.8 is relatively self-consistent, since the short-distance 
contributions begin to dominate. The above percentages of perturbative 
contributions are smaller than those in the heavy meson decay B — > ttIu, 
which are roughly 60-70% near the high end of the pion energy spectrum 
[|l|]. This is expected, because there are more partons to share the baryon 
momentum, such that the partons are softer, and Sudakov suppression is 
weaker. 

Compared to the analysis in [|13|], we have taken into account the resum- 
mation effect from the A^ baryon and the contributions to the hard subampli- 
tude from Figs. l(g)-l(n). Recalculating Li for (3 = 1.0 and p = 1.0 without 
the Sudakov factor associated with the A& baryon, we observe a 37% decrease. 
Note that the inclusion of the Sudakov factor from the heavy baryon causes 
an enhancement, instead of a suppression. It is not difficult to understand 
this result. Since in Eqs. ( |26| ) and fl2T| ) are not positive definite in the 
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whole kinematic region, it is possible that the modulation of this Sudakov 
factor suppresses the negative portion of the integrand in Eq. (|23|) more than 
the positive portion. 

The more essential improvement is the consideration of Figs. 1(g)- l(n) 
in the evaluation of the hard subamplitude. It was argued that the con- 
tributions from these diagrams, with at least one of the exchanged gluons 
attaching the b quark, are suppressed by the power 1/M^ b , and thus negligi- 



ble jnj. However, this argument is correct only when the exchanged gluons 
are soft, such that the b quark propagators in Figs. l(g)-l(n) provide a power 
suppression by 1/M\ b . Before investigating the validity of this argument, it 
has contradicted the assumption of hard gluon exchangs, based on which the 
PQCD formalism was developed. Our analysis indicates that with the help 
of Sudakov suppression, the exchanged gluons are indeed hard enough, and 
that Figs. l(g)-l(n) lead to contributions of the same order as those from 



Figs l(a)-l(f). If considering only the first six diagrams in Fig. 1 as in [13 
the result of L\ for (3 = 1.0 and p = 1.0 reduces by a factor 3. That is, the 



diagrams ignored in [ 1 3 1 are in fact even more important. 



At last, we present the variations of L\ and of the differential decay rate 

% = I^M 3 "^^ (36) 

for massless leptons with the proton energy fraction p. Employing the Fermi 
coupling constant Gf = 1.16639 x 10 -5 GeV -2 and the CKM matrix element 
V u b = 0.003, results for (3 = 1.0, 1.5, and 2.0 are displayed in Figs. 5 and 
6, respectively. The curves rise quickly as p < 0.6, implying the dominance 
of nonperturbative contributions. Our predictions become insensitive to the 
parameter (5 of the baryon wave function as p > 0.8. Therefore, we 
conclude that the values of L\ and of dT /dp for p = 1.0 are about L\ = 
2.3 x 10~ 3 and dT/dp = 9.4 x 10~ 21 GeV. 



VI. CONCLUSION 

In this paper we have studied the semileptonic decay — > pW as an 
initial application of the PQCD factorization theorem to heavy baryon de- 
cays. It is found that the perturbative analysis is reliable near the high end 
(p > 0.8) of the proton energy spectrum. Compared to the previous study 
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, we have taken into account the Sudakov resummation for a heavy-light 
system (the A& baryon) and a complete hard subamplitude. The quark-level 
decay diagrams neglected in |T3| turn out to be very important. Our approach 
will be applied to the semileptonic decay A& — > A C W, for which experimental 
data are available. We shall try to determine the parameter (3 of the A& 
baryon wave function from this decay, and then extend the PQCD formalism 
to nonleptonic A^ baryon decays, which are more challenging. These subjects 
will be published elsewhere. 
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Table 1. The expressions of \I/ pa'p'-y' H a ' /3 ' y ' al3 ' 7 ^ A b a/3-y from the diagrams 
in Fig. 1 in terms of 

Va = Xix' iP Ml b + (k iT - Kt? , 

For each diagram, (a)-(n), row (1) is the original expression, row (2) is the 
series of variable changes necessary for bringing row (1) into the desired form, 
and row (3) is the expression after variable changes. 



Diagram ^ P H^ Ab /(47i 2 a 2 s f P f A j27) 



1^213^123 +2Tl230123j ~ 

(l-x 1 ){l-x , 1 )V 11 V 33 

[0312<ftl32+2Tl32 0132] 

(l-xi)(l-xi)X»nX»22 






2Tl23 0*123 

(l-xi)(l-x' 2 )V 22 V 33 

l->3,2->l,3->2 

2T 3 i2^3l2 

(l-x 3 ){l-x')V 11 V 22 

— '02139123 

(l-x 1 )(l-x' 3 )V 22 V 33 

l->3,2->l,3->2 

— ^1320312 

(l-X 3 )(l-x' 2 )X»iiX»22 





["02130123 +271230123] 
(l-Xl)(l-xi)X»uX»22 

None 

[•02130123 +271230123] 
(l-Xl)(l-x' 1 )X'iiX'22 
2x 3 (p~l)?123 0123 
(l-X 2 )p[(l-X3)x 2 p+X 3 ]X>22l>33 

l->3,2->l,3->2 

2x 2 (p— 1)T312</>312 

(1-3^ )p[{l-X 2 )x\ p+X 2 ]PilX>22 



(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



1) 
2) 
3) 

1) 
2) 
3) 
1) 
2) 
3) 

1) 
2) 
3) 

1) 
2) 
3) 
1) 
2) 
3) 

1) 
2) 
3) 
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Diagram ^pH^ Ab /(A7r 2 a 2 J P f Ab /27) 

OO (i) o 

(2) 

(3) 



(i) (1) 

(2) None 



(3) 
(j) (1) 



(3) 

00 (i) 



(3) 
0) (1) 



(3) 
(n) (1) 



(3) 



[V'2130123 +2Ti230123] 

(i-x' 1 y P v 11 v 2 2 



x'l ["0213 0123 +2Ti23 0123] 
(l-xi) 2 pX>iiX>22 
[V>213 0123 +2Tl23 0123] 
(l-Xl)(l-xi)X>iiX>22 

(2) None 

[•02130123+27123 0123] 
(l-Xl)(l-x' 1 )X>nX>22 
[(l-X3-a / 1 a2)p-X2]-02130123 
(l-xi)(l-X 2 )p(x^/3-l)X»22X>33 
2[(l-2x2+x 3 -x 3 p)p+X2]Ti230i23 
(1-X' 1 )(1-X 2 )/>(X^-1)X>22X>33 

(2) l->3,2->l,3->2 

[(l-X2~X33:i)p-Xl]V>l320312 
(l-X^)(l-Xl)p(x' 2 /3-l)X>llX>22 

2[(l-2xi+x 2 -x 2 p)p+xi]T 3 i 2 <fel2 

(l-X^)(l-Xl)p(x' 2 /9-l)X»llX>22 
[(l-2x 3 +X^-X 2 p)p+X 3 ]^2130123 



+ 



+ 



(l-x')(l-X 3 )p(x 2 p-l)X>22X>33 
2[(l-x 3 -a; 2 -(l-3;3)(l-x / 1 )p)p-3:3]Ti230i23 
~T (l-x' 1 )(l-X3)p(x' 2 p-l)X» 22 X>33 



(2) l->3,2->l,3->2 
(3) 



[(l-2x 2 +x' 1 -X / 1 p)p+X2]V'13203 
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(l-x')(l-x 2 )p(x' 1 p-l)X'iiX> 22 

I 2[(l-X2-X 1 -(l-X 2 )(l-X 3 )p)p-X2]T 3 i20312 
(l-X 3 )(l-X2)p(x' 1 p-l)D 1 iI?22 
(' rn \ /1\ [ - 02130123+2Ti230l23] 
^ ^ (l-Xl)(l-xi)X>iiX>33 

(2) 2^3 



[^312 0132+27132 0132] 
(l-Xl)(l-x' 1 )X>nX>22 
x[ [V'2130123 +2Tl230123] 

(2) 2^3 

ft'l ["03120132+2Ti 3 20132] 
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Table 2. The constants and the Appel polynomials involved in the 
proton wave function if) in Eq. (|30"D |22| . 



3 % Nj b j A A x i) 

1.00 1 1 

1 0.310 21/2 20/9 zi-zg 

2 -0.370 7/2 24/9 2 - 3(xi + x 3 ) 

3 0.630 63/10 32/9 2 - 7(zi + x 3 ) + 8(x? + s§) + ± Xl x 3 

4 0.00333 567/2 40/9 x x - x 3 - (4/3) {x\ - x\) 

5 0.0632 81/5 42/9 2 - 7( Xl + x 3 ) + U Xl x 3 
+(14/3)(xf + x§) 
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FIGURE CAPTIONS 

FIG. 1 Lowest order diagrams for the A b — > pW decay. 
FIG. 2 Radiative corrections to Fig. 1. 

FIG. 3 Dependence of L\ on b c for (3 — 1.0. b c is in unit of 1/Aqcd- 
FIG. 4 Dependence of L x on b c for (3 = 2.0. b c is in unit of 1/Aq C d- 
FIG. 5 Dependence of L\ on p. 

FIG. 6 Dependence of dT /dp on p in unit of 1CT 18 GeV. 
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